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ABSTRACT 

In addition to the dominant oscillatory gravitational wave signals produced during binary inspirals, a non- 
oscillatory component arises from the nonlinear "memory" effect, sourced by the emitted gravitational radia- 
tion. The memory grows significantly during the late inspiral and merger, modifying the signal by an almost 
step-function profile, and making it difficult to model by approximate methods. We use numerical evolutions 
of binary black holes to evaluate the nonlinear memory during late-inspiral, merger and ringdown. We identify 
two main components of the signal: the monotonically growing portion corresponding to the memory, and an 
oscillatory part which sets in roughly at the time of merger and is due to the black hole ringdown. Counter- 
intuitively, the ringdown is most prominent for models with the lowest total spin. Thus, the case of maximally 
spinning black holes anti-aligned to the orbital angular momentum exhibits the highest signal-to-noise (SNR) 
for interferometric detectors. The largest memory offset, however, occurs for highly spinning black holes, 
with an estimated value of /ijo' — 0.24 in the maximally spinning case. These results are central to determin- 
ing the detectability of nonlinear memory through gravitational wave interferometers and pulsar timing array 
measurements. 

Subject headings: black hole physics — gravitational waves 



1. INTRODUCTION 

Accelerating massive bodies generate gravitational waves 
(GWs), a potentially rich source of astrophysical information 
which a number of large experiments have been designed to 
mine over the next decade. The principle sources for the inter- 
ferometric detectors (LIGO, Virgo, LISA) are orbiting bodies, 
for which the dominant signal is oscillatory. However, it has 
long been known that GWs will also contain non-oscillatory 
features, resulting from a net chang e in time derivatives of 
the m ultiple moment s of the system ( jZel'Dovich & Polnarev| 
1974 jBraginskii & Thorne 1987| l. A nonlinear contribu- 



tion also results from the interaction of the waves with them- 



selves ( Payne 1983 Christodoulou 1991[ Blanchet & Damour 
[1992). Evaluating these modes involves an integral in time 
requiring knowledge of the entire past history of the space- 
time. As such, they have been given the name "memory" or 
"hereditary" components of the GW. 

The nonlinear memory is sourced by the emitted GWs. 
For a relativistic binary inspiral, the typical profile is of a 
slow growth over time which sees a rapid increase during 
the late inspiral and merger, and reaches a constant value 
as the merger remnant rings down and ceases to emit. The 
non-oscillatory nature of the memory suggests that it will not 
be a dominant feature in interferometric detectors, for which 
drifts in the background metric are factored out. However, 
the step induced during merger may be observable in pulsar 
timin g measurements, as has been noted in a recent set of pa- 
pers tPshirkov et al.|2009Hvan Haasteren & Levin|2009[ [Setol 



[2009 



Post-Newtonian (PN) calculations provide an accurate esti- 
mate of the m emory to within a few orbits of the merger ( KeT>] 
nefick 1994| l, however fail to model the important merger 



'denis.pollney@uib.es 
^reisswigigtapir.caltec n.edul 



phase where the e ffect is largest. In a series of papers, jFavata 
( 2009a|c|b 2010) has included an estimate of the merger con- 
tribution for binary black holes (BHs) by using an effective- 
one-body (EOB) model which was tuned using the results of 
numerical simulations. Recent progress in numerical models 
of BH spacetimes provides the op portunity to mea sure the ef- 
fect directly, though as outlined in |Favata ( |2009c| l, due to the 
low amplitude of the memory, and systematic error in stan- 
dard techniques of wave measurement in numerical relativity, 
the problem is a challenging one. 

In t his paper, w e apply a newly developed evolution 
code ( [Pollney et al . 2009b a) and wave extraction tech- 
niques ^Reisswig eral.„2009a, ,2010 ) to carry out fully rela- 
tivistic binary BH simulations to directly measure the nonlin- 
ear GW memory through inspiral, merger and ringdown in a 
gauge-invariant and mathematically unambiguous notion. We 
concern ourselves with equal-mass, spinning binaries. For 
these models, the memory is contained in the -2Ym spin- 
weighted spherical harmonic components of the GW strain. 
These m odes, previously stud ied in the context of head- 
on collisions ( [Anninos et al.|1995] l, exhibit two dominant ef- 
fects. The first is a non-oscillatory term which rises as the 
inspiral progresses, and is associated with the memory. Dur- 
ing the merger, an oscillatory signal is superposed onto these 
modes, induced by the ringdown of the BH remnant. After 
the ringdown has subsided, the modes are offset from their 
original value. We find that for the dominant {£, m) = (2,0) 
memory mode, this offset is largest in the case of the merger 
of aligned, maximally spinning BHs, as might be expected 
given that these are the strongest gravitational emitters ( jReis^ 
[swig et aL][2009b) l. Interestingly, we observe the strongest 
oscillatory ringdown signal in the case of lowest total spin 
(maximally spinning BHs anti-aligned with the orbital angu- 
lar momentum). As such, and as opposed to the observability 
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of the dominant {£, m) = (2, 2) mode, these models produce 
a more visible (£, m) = (2, 0) mode to interferometric detec- 
tors. We conclude by estimating the strain offset which would 
be visible to a pulsar timing array. 

2. NUMERICAL METHODS 

We integrate the BSSN OK formulation of the vacuum Ein- 
stein equations (see, e.g., ' Alcubierre 2008 ) numerically, us- 
ing the Llama evolutio n code (Pollney et al. 2009a b) within 
the Cactus framework ( |Gooda le et al. 2003 ). The evolution 
equations are discretized via finite differences in space on a 
grid composed of multiple patches with locally adapted co- 
ordinates in the wave zone. Adaptive mesh-refinement is im- 
plemented via the Carpet grid driver (|Schnette r et al.|[2004 



on the central patch surrounding the BHs. Importantly for 
the measurement of small non-oscillatory features such as the 
memory, the artificial grid outer boundary is causally discon- 
nected from the wave measurements. 

We determine initial dat a for binary BHs by the confor- 
mal puncture data method ("Brandt & Briigmann 1997), and 
evolve them using the "moving puncture" approach ( Baker] 
et al.| (2006, Campanelli et al. 2006) whereby the choice of 



gauge ( [Alcubierre et al.,2003^ prevents the sp acetime slicing 
from encountering the curvature singularity ( |Hannam et al. 
[2007] l. 

Crucial for determining the memory, we compute the GW 
signal of the source using the method of Cauchy-characteristic 
extraction (CCE), which computes unambiguous and coordi- 
nate invariant signals at future null infinity, J'^ , correspond- 
ing to a detector far removed from the source (Bishop et al. 



1999; Babiu c et aL|[200 5; Winicour 2005; Reisswig et al. 



5D09a, 2010). By this technique, metric data is collected on a 
world-tube at finite radius, and a formulation of the full Ein- 
stein equations along null hypersurfaces is used to transport 
the signal to . 

We measure two quantities that encode the gravitational 
signal at . The first is the Newman-Penrose scalar ^4, 
the Weyl curvature component with slowest falloff in asymp- 
totically flat spacetimes (Newman & Penrose 19 62]l. Alter- 



natively, we measure the Bondi "news" (Bondi et al. 1962 
Sachs|1962| l, defined by 



N = -Aa, 



(1) 



where A — V^V a is a derivative operator defined on an out- 
going null geod esic, and a is th e Newman-Penrose shear 
scalar ( Newman & Penrose||1962] ). These scalars, t/;4 and 
M, can be evaluated directly from the local curvature on the 
sphere at J" "'"a nd in the naturaUy defined gauge of an asymp- 
totic observer ( |Bondi et aL]|1962| |Sachs||1962| l. To calculate 
the GW strain h which is observed by a detector, these quan- 
tities need to be numerically integrated in time: 



h = — ihy 



dt'N ^ 



dt' / dt"'(/'4- (2) 



We expand h in terms of a basis of spin-2 spherical harmon- 
ics -2Yem- The dominant non-oscillatory hereditary compo- 
nent is contained in the {£, m) — (2, 0) mode. 



(3) 



r/i2o = / (iri(_2>2o)/j, 
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Fig. 1. — Upper panel: the {I, m) = (2, 0) mode of the news, Af. Lower 
panel: the strain, h (integrated from A^). The dotted line in the inset of the 
upper panel is the result of a fit to the ringdown QNM, Eq. |5] while in the 
lower plot it denotes the 3PN memory deiived in liFavata (i20b9c), used to 
determine the integration constant. 

Since the simulations are necessarily of finite length, the 
time integration of M, Eq. (|2]), leaves a constant to be de- 
termined, corresponding to the initial value of the signal at 
the start time of the simulation. For purely oscillatory signals 
(e.g., /122) we can fit this by adjusting the post ringdown value 
to zero, or by some averaging procedure over multiple inspiral 
wavelengths. This is not possible in the case of the monoton- 
ically growing memory. Instead, we approximate an initial 
value of the ft,20 mode by matching to a PN calculation, using 
the 3PN expansion for derived inlFavata (2009c ). For 



where J7 is the sphere at J'"'" and r is the distance to the source. 



spinning binaries we incorpor ate 2.5PN spm contributions to 
the frequency evolution (Blanche! et al. 2006] l. The initial off- 
set of the numerical waveforms is performed by shifting its 
amplitude so that it fits against the PN models, over an inter- 
val from t = -300A/ to t = -200A/. For instance, the PN 
solution for /120 in the non-spinning case is plotted as a dot- 
ted line in the lower panel of Figure [T] In Table [T] we report 
both the overall amplitude offset of the ringdown signal in the 
/120 mode matched to the PN inspiral, as well as the offset 
from the point t — — 300M which is computed entirely from 
numerical data. 

3. CHARACTERISTICS OF THE MEMORY MODES 

In Figure [1] we show the {i,m) — (2,0) mode during 
late inspiral and merger for the case of an equal-mass, non- 
spinning b i nary us ing data from the model p resented in [Poll 
Iney et al. ( 2009b 1 and evaluated at in Reisswig et al. 
(2009aj). ITie upper panel shows the development of the 
{i,m) = (2,0) component of the news, Af2o, whereas the 
lower shows its integral, /i20- The time coordinate places 
t — OA/ at the position of the amplitude peak of the dominant 
(2, 2) mode (not shown), corresponding roughly to the time 
of merger. Two main features are present. The first is a slow 
monotonic growth of the A/20 mode during the inspiral and 
plunge phase. This corresponds to the non-oscillatory mem- 
ory contribution. At t ~ — 20Af , there is a local maximum in 
JV20, and shortly afterward a prominent oscillatory signal sets 
in. This signal represents the binary merger and subsequent 
BH ringdown. It is not strictly an aspect of the GW "mem- 
ory", as it is not dependent on the entire past history of the 
spacetime, however it is superposed on the growing memory 
signal. 
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The lower plot shows the strain, determined via Eq. l2] The 
superposition of the memory and ringdown signals results in a 
notable kink in the waveform near t = QM, and exponentially 
decaying oscillation before reaching the steady-state ampli- 
tude of 

/i'2°o = 0.097±2 X 10"^. 



(4) 



We measure approximately 4-th order convergence in A/20 
when comparing three different grid resolutions for the same 
model, and we estimate an error on the order of 1 x lO^"* in 
Af2o- The much larger error of 2 x 10^"^ in is an upper 
limit of the estimate on the total error composed of integration 
error (2]) and error in the PN fit for all of the binary models 
considered in this paper 

Perturbative results provide a pred iction for the obs erved 
quasi-normal modes (QNM) of A/'2o ( |Berti et al.]|2009| l. We 
perform a least-squares fit of the numerical results to a func- 
tion of the form 



f{t) — Aexp(— t/r) cos(tj2oi + > 



(5) 



over the ring-down portion, t E [30, 80], with fitting parame- 
ters, A, and W2o^ 4>^ ™d arrive at a frequency 



MfW2o = 0.3940 ± 2 X lO"" , 



(6) 



using the measured mass of the remnant of Mf = .951764 ± 
20 X 10"^ ( jPoUney et al |2009b| ). [Beiti et al.| ( |2009 l provide a 
tabulated value of Muj^^^- = 0.393245 for the (2, 0) mode of 
a BH with the measured dimensionless spin Of = 0.686923 ± 
1 X 10^^. Thus, we find a difference of 



lAffCJao - Moj^qI « 8 X 10 



(7) 



representing an error of less than 0.2%. We list the final spin 
and the associated QNM frequencies of the binary systems for 
all simulations performed in Table [T] 

Finally we have also examined m = modes for £ ~ 4 
and £ = 6. For this model, the overall amplitudes of the fi- 
nal memory are h^^^ = 1.7 x 10"^ and /i^g^* w 5.0 x 10""*. 
However, these modes are difficult to distinguish from the nu- 
merical error, and we conclude that the higher order modes 
contribute less than 2% to the total memory. 

4. MEMORY FROM SPINNING BINARIES 

We have performed simulations using the code infrastruc- 
ture described in Section|2]in the two-dimensional parameter 
space of non-precessing equal-mass binaries with spin that 
has also been studied in Reisswig et al. (2009b); Rezzolla 



et al. (2008 1. We have focused on initial data configurations 
along the two main axes in the space of individual dimen- 
sionless BH spins (01,02) aligned with the orbital angular 
momentum. The first set of models have equal but opposite 
spins, fli = —02; The second have identical spins, ai — a2- 
The model details are listed in Table [T| 

Figure l2]plots the evolution of the /i2o modes for the two se- 
quences. The upper panel shows models for which the spins 
of the individual BHs are anti-aligned, ai = —02, so that 
the spacetimes have the same total angular momentum. The 
nearly identical lines are consistent with previous observa- 
tions that the {£, m) = (2, 2) modes also do not vary appre- 
ciably alon g this direc tion of parameter space ( Vaishnav et al. 
|2007, Reis swig et al . 2009b), as well as the radiated ener- 
gies (Reiss wig et al.|[2009bj ), at least not to within the error- 
bars as stated in the latter reference. 
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Fig. 2. — Upper panel: the (I, m) = (2, 0) modes are similar for models 
with anti-aligned spins, ai = — a2.n Lower panel: the memory offset grows 
with the total spin, however the ringdown is most prominent for the anti- 
aligned spins, (ai, 02) < (0, 0). 

The lower plot shows the /i2o modes for configurations with 
equal and aligned spins, ai = 02. The amplitude of the mem- 
ory offset increases with higher total spin. This is expected as 
the memory is sourced by the emitted GWs w hich are known 
to be more energetic in the higher spin models ( Reisswig et al. 
2009b; Loustoetal.|2009 ). The total offsets, yq^, are listed in 
Table [T] 

We find that the amplitude of the ringdown is largest for 
models with lowest total spin, ai = a2 ~ —0.8. In these 
cases, spin-orbit coupling increases the final angle of impact 
so that the merger is more nearly head-on. The ringdown is 
suppressed in the high spin ai = 02 = +0.8 model. Fits 
to the QNM frequencies are listed in Table [T] and show good 
agreement with perturbative results. However, for simulations 
with ai = 02 > the weakness of the ringdown contributes 
to larger variation in the estimates, particularly if the fitting 
interval is varied, though still within 5% of the analytic values. 

The memory is well estimated by a quartic polynomial in 
the total spin a = (ai + a2)/2, 



'20 



(a)r : 



0.0562 a + 0.0340 + 0.0296 



0.0206 a'' 



(8) 



determined by a least-squares fit to the measured data. Fig- 
ure [sjplots /ijo' as a function of spin for the aligned models, 
fli = a2, as well as the leading order multip ole moment esti - 
mate resulting from the radiated energy (see |Favata (2009a), 
Eq. (5)) 



h^^ ^ ^^^^"J sin^ e{l7 + cos'' 



with 



3847ri? 

^mem 



16^ ( AE,,i 



(9) 



(10) 



The measured radiated energy, A_Er ad, can also be estimated 
using the quadratic fit developed in [Reisswig et al. | ( |2009b[ ), 
verified here with the plotted data points corresponding to the 
new and more accurate simulations (see Table [TJ. The radi- 
ated energy provides a reasonably good estimate of the mem- 
ory offset, though with a slight under-prediction. The mea- 
sured memory values also appear to have a somewhat stronger 
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Fig. 3. — The total memory, h^g, given by the offset of the /120 after ring- 
down. The solid line is a quartic fit through the measured data (+) with 
integration constants determined by PN. This is compared with the estimate 
arising from the total radiated energy, via Eq.|9]( X , dashed line). 

dependence on the total spin. Extrapolating to the extremal 
ai = 02 = 1 case, we estimate a maximum /iJfo — 0.24 in 
the case of aligned spins. Fits to AE^.^^ for generic spin con- 
figurations (Lousto et al. 2009) suggest that this will indeed 
be the maximum value for arbitrary spins. Additional simu- 
lation in the high-spin regime would be required to establish 
the accuracy of the extrapolation. 

5. DETECTABILITY OF THE MEMORY MODES 

The m = modes are small (on the order of 10%) com- 
pared to the dominant £ = m components of the GW sig- 
nal. However, they exhibit both non-oscillatory growth and 
ringdown features, which present some interesting aspects for 
detection and identification of these modes. We discuss the 
prospects for observation in both interferometers, and pulsar 
timing arrays, in the following sections. 

5.1. Interferometers 
The signal-to-noise ratio (SNR) within a given detector, p, 



in GW searches based on matched filtering is given by Flana- 
gan & Hughe"s|([l998|): 



„2 _ 



= 4 



matched 



\Hf)l 

Shif) 



df, 



(11) 



where Sh{f) is the sensitivity curve of the instrument. Re- 
stricting the integral to the {£, m) = (2, 0) mode, we can de- 
termine its contribution to t he o verall SNR, testing against the 
propos ed advanced-LIGO ([Shoemaker 2010) and LISA noise 
curves (LISA-wiki 2010). 

Table [T] lists the maximum SNR attained over a range of 
masses for a given model. For advanced-LIGO, the optimal 
total mass is Mbh — 29OM0 (independent of the model and 
invariant under the distance), i.e. a 145M0 + 145M0 bi- 
nary, for which we list the SNR at a reference distance of 
SOOMpc (the SNR scales linearly with the distance). The 
LISA results refer to an optimal (redshifted) total mass of 
AiflH = 5.35 X IQ^Mq at a luminosity distance d — 15.8Gpc 
corresponding to a redshift z = 2, wher e models suggest a t 
least one expected merger event per year ( [Sesana et al. 2005| l. 
By assuming a minimum SNR of p = 8 for detection, it will 
be possible to observe the (£, m) = (2, 0) mode in LISA out 
to the given distance. They are unlikely to be visible with 
significant event rates within the planned advanced ground- 
based detectors, though future generation experiments may 
have prospects. 

For both detectors, and as opposed to the dominant 
{£, m) = (2, 2) mode ( [Reisswig et al.||2009b| ), we note that 
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Fig. 4. — Fourier-representation of (£, m) = (2, 0) for different equal and 
aligned spin models at the optimal total mass M = 290Mq (145M0 + 
145A/0 binary) at a distance d = SOOMpc and at an angle 9 = tt/2. 
Models with higher total spin (e.g. ai = 02 = +0.8) have a weaker ring- 
down oscillation, which results in a higher overall slope and thus smaller 
amplitude at high frequencies as opposed to models with lower total spin 
(e.g. ai = 02 = -0.8). 

the highest SNR occurs for the cases with spins anti-aligned 
to the orbital angular momentum which is due to their stronger 
ringdown. In the aligned cases, the weaker ringdown causes a 
much flatter high-frequency fall-off (Figure]?]), and this effect 
overrides the non-oscillatory step within these detectors. 

An important aspect of the matched-filtering algorithm is 
the construction of templates against which signals are com- 
pared. We examine the influence of neglecting the (2. 0) 
contribution by computing the mismatch between a signal, 
h{6 — 7r/3,(/) ~ 7i'/2), containing all of the (2,to) modes, 
with one in which the (2, 0) mode is left out. For advanced- 
LIGO, the mismatch for the model ai — a2 = +0.8 is largest 
in the mass range 25OM0 — 4OOA/0 , but negligible, with a 
magnitude of 10^^^. For the ai = 02 = —0.8 model, the mis- 
match grows to be on the order of 10^^. Thus, although the 
ft,20 mode induces a notable offset in the wave signal when 
viewed in the time domain, it makes little difference to the 
detection algorithms. 

5.2. Pulsar timing arrays 

An alternative proposal for measuring GWs may be more 
sensitive to the non-oscillatory step-function nature of the 
memory. Precise, distant clocks will experience a residual in 
time-of-arri val measurements which can be associated with 
GWs ( Estabrook & Wahlquist|1975 [lSazhin'1978','Detweiler 
1979 ) . Several such projects are incl evelopment (Hobbs et aT 
|20D9||Manchester|2008||Jenet et al.|2009nLazio|2009f While 
their principal aim is to examine the stochastic background, 
recently a number of papers have pointed out the potential 
for observing n on-oscillatory step function burst signals, such 
as the memory ( Pshirkov et al.|20()9{ van Haasteren & Levin 



]2009l]Setol2009l l 

Primary candidates for detection are supermassive BH bi- 
nary mergers. While there is a great deal of ongoing work 
understanding the nature and evolution of such systems, if 
there is a tendency for the bodies to align, then the results 



obtained here are directly applicable. For instance, [Dotti et aL 
( 2009[) indicate that in gas-rich mergers, the spins are aligned 



on a short timescale, with spin magnitudes in the range 0.6 
to 0.8. Assuming a mi = m2 = lO*Af0 binary source at a 
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TABLE 1 

Parameters and measurements of physical quantities for the equal-mass aligned-spin BBH models. From left to right, we list: 

INITIAL DIMENSIONLESS SPINS ai AND a2, THE TOTAL MEMORY /l^'g* ACCUMULATED FROM T > — SOOAf PRIOR TO MERGER AS COMPUTED FROM NR 
simulations, and THAT ACCUMULATED FROM T > -OO BY MATCHING PN TO NUMERICAL DATA, THE MAXIMUM SNR p OF THE (2, 0) COMPONENT 
ATTAINED OVER A RANGE OF MASSES (Mopt = 290Mq FOR ADLIGO, Afopt = 5.35 X 10^ FOR LISA) AT LUMINOSITY DISTANCES d = 300Mpc FOR 
ADLIGO and d = 15.8Gpc (z = 2) for LISA, where the bracketed numbers report the sky-averaged SNR, THE TOTAL RADIATED MASS 
Brad, INCLUDING CONTRIBUTIONS FROM PN, AS WELL AS THE FINAL DIMENSIONLESS SPIN afin OF TH E REMNANT, AND FINA LLY, THE MEASURED 
QNM FREQUENCIES McJ20 AND THE ERROR TO PERTURBATIVE CALCULATIONS 0F |BERTI ET AL.| ( [2009l l. 



ai 02 


'^20*lTg[-300,1001 ^^20 


(mem) (mem) 
PadLIGO PUSA 


Srad (%) 


"fin 


Mu)20 |Ma;Jit' - A/lJ2o| 


-0.2 0.2 
-0.4 0.4 
-0.6 0.6 
-0.8 0.8 


0.073 0.097 
0.073 0.097 
0.073 0.097 
0.072 0.095 


5 (4) 62 (51) 
5 (4) 63 (52) 
5 (4) 65 (54) 
5 (4) 78 (64) 


4.84 
4.85 
4.91 
4.94 


0.6866 
0.6865 
0.6844 
0.6832 


0.3932 5 X 10-4 
0.3932 5 X 10-4 
0.3931 3 X 10-4 
0.3978 48 X 10-4 


-0.8 -0.8 
-0.6 -0.6 
-0.4 -0.4 
-0.2 -0.2 
0.0 0.0 
0.2 0.2 
0.4 0.4 
0.6 0.6 
0.8 0.8 


0.044 0.067 
0.049 0.072 
0.056 0.078 
0.064 0.084 
0.074 0.097 
0.085 0.110 
0.102 0.126 
0.126 0.153 
0.162 0.188 


11(9) 147(121) 
9 (7) 128 (106) 
7 (6) 103 (85) 
6 (5) 81 (67) 
5 (4) 63 (52) 
4 (3) 47 (39) 
3 (2) 36 (30) 
2 (2) 28 (23) 
1(1) 21(17) 


3.42 
3.66 
3.97 
4.35 
4.83 
5.44 
6.22 
7.33 
8.93 


0.4247 
0.4919 
0.5604 
0.6246 
0.6869 
0.7469 
0.8041 
0.8575 
0.9039 


0.3824 19 X 10-4 
0.3824 6 X 10-4 
0.3850 11 X 10-4 
0.3905 11 X 10-4 
0.3940 8 X 10-4 
0.3987 12 X 10-4 
0.4010 13 X 10-4 
0.4080 6 X 10-4 
0.4035 89 X 10-4 



distance of IGpc which has aligned spins with a magnitude 
of ai — a2 ~ 0.7, the results of Section]?] lead to an angle- 
averaged memory offset of: 

(to%)(1")- 

(The corresponding constants for zero and maximally spin- 
ning bodies ar e 0.9 x IQ-^^ and 2.0 x 10 -^^, respectively.) 
The results o f [Pshirkov et al.| ( |2009| ) and [van Haasteren &| 
|Levin| ( |2009 ^) indicate that a burst will be observable over a 
10-year observation period at amplitude /i ~ 2 x IQ-^^, sug- 
gesting that the merger of a pair of 6 x 10^ with ai = 02 = 0.7 



would be observable at IGpc. 
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